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We present a new version of the CompHEP program (version 4.4). We describe shortly new issues implemented 
in this version, namely, simplification of quark flavor combinatorics for the evaluation of hadronic processes, Les 
Houches Accord based CompHEP-PYTHIA interface, processing the color configurations of events, implemen- 
tation of MSSM, symbolical and numerical batch modes, etc. We discuss how the CompHEP program is used 
for preparing event generators for various physical processes. We mention a few concrete physics examples for 
CompHEP based generators prepared for the LHC and Tevatron. 



1. INTRODUCTION 

The increase of collider energies requires effi- 
cient and precise calculation of processes with 
multi-particle final states. The NLO, NNLO or 
resummation corrections to the leading order re- 
sults should be incorporated when available. 

At LEPl the most interesting physical chan- 
nels were 2-fermion production processes at the Z 
pole, while at LEP2 mainly 4-fermion processes 
with some photons or gluons radiated off have 
been considered. At the Tevatron, LHC and Next 
Linear Colliders a majority of interesting channels 
include 5,6 and 8 or even more fermions in the fi- 
nal states with additional photons or gluons. For 
example, the top quark pair production leads to 
6 fermions in the final state, single top produc- 
tion in various modes gives 4, 5 and 6 fermions, 
the channel till produces 8 fermions in the final 
state, etc. Calculations of event characteristics 
with multi-particle final states are needed to be 
done at the complete matrix element level in or- 
der to include correctly the spin correlations, non- 
trivial kinematics and in some cases finite width 
effects in a gauge invariant manner. 

In practical computations, especially for 
hadron colliders, there are many contributing 
Feynman diagrams even at the tree level and 



a large number of subprocesses leading to the 
same multi-particle final state. To perform such 
computations a number of automatic programs 
have been created - CompHEP [Tj, GRACE H], 
MadCraph (31, AlpGen g], O' Mega with 
WHIZARD 6 , Amegic 0, HELAC-PHEGAS 
and other. 

However, in order to get reliable predictions 
and reduce a dependence on scales involved, like 
QCD normalization and factorization scales, one 
should compute and simulate not only processes 
with many final particles but also include various 
corrections and perform resummations of large 
logs. Recently a few calculations of NLO cor- 
rections for 2 — > 3 processes have been performed 

Precise calculations including the NLO correc- 
tions to the rate and distributions are however not 
sufficient for experimental studies insofar as one 
should also provide a corresponding NLO event 
generators. Since the real final states are not par- 
tons but hadrons (jets) these generators should be 
matched in some way to the parton shower gen- 
erators, such as PYTHIA [TS| and HERWIG [Tl] . 
which simulate a realistic picture of hadroniza- 
tion. In fact the problem of matching of all the 
needed ingredients, multi-particle matrix element 
calculations with all spin correlations, NLO cor- 
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rections to, at least, basic production processes, 
hadronization taking into account initial and final 
state radiation is highly nontrivial. As a result 
one likes to have an event generator which in- 
cludes all these parts without "double counting" 
of contributions and with smooth distributions in 
the whole range of the phase space. Although this 
problem is far from a final solution, recently in- 
teresting and promising approaches for different 
parts of the above mentioned matching have been 
proposed tL5 16 17.18..19-20, . 

In this paper we discuss mainly a new aspects of 
the CompHEP latest version 4.4 and present some 
recent CompHEP applications to physics studies 
at colliders. We discuss shortly the SingleTop 
generator, created with the help of CompHEP to 
simulate the electroweak top quark production at 
the Tevatron and LHC, where proper matching 
of different processes with main contributions to 
different phase space regions is performed in the 
final event flow. 

CompHEP is based on quantum theory of 
gauge fields and includes the Standard Model 
Lagrangian in the unitary and t'Hooft-Feynman 
gauges as well as several other MSSM based 
models. CompHEP is able to compute basically 
the LO cross sections and distributions but with 
many particles (up to 4-6) in the final state tak- 
ing into account, if necessary, all the QCD and 
EW diagrams, masses of fermions and bosons and 
widths of unstable particles. Processes computed 
by means of CompHEP could be used as a " new" 
external processes for generators like PYTHIA. 

The CompHEP project started in 1989 at SINP 
MSU. During the 90's the package was developed, 
and now it is a powerful tool for automatic com- 
putations of collision processes. The CompHEP 
program has been used in the past for many stud- 
ies (see |2I] and the citation to CompHEP 1, for 
more complete information) as shown schemati- 
cally in the Fig^ 

2. COMPHEP 4.4 

In the latest CompHEP version 4.4 
(Ihttp: / / theory. sinp . msu . ru / comphep ) there are 
several new improvements incorporated which 
make the program significantly more efficient to 
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Figure 1. Very incomplete list of processes com- 
puted by means of CompHEP in the past 

use for event simulation, especially for the case 
of hadronic collisions. 

2.1. QUARK FLAVOR COMBINA- 
TORICS 

A serious computational problem for the 
hadronic collisions is a large number of partonic 
subprocesses with many contributing Feynman 
diagrams in each of them. This takes place, 
in particular, for the reason of numerous possi- 
ble combinations of initial and final quarks with 
different flavors, giving different subprocesses. 
Many additional diagrams appear in each of these 
subprocess due to CKM quark mixing. In fact 
many of the subprocesses and many diagrams 
have very similar topological and Lorentz struc- 
ture, so one would like to compute all the similar 
pieces only once. 

In the paper 153' the method for a simplifica- 
tion of quarks combinatorics has been proposed. 
The method has been implemented to the Com- 
pHEP 4.4 and tested on many examples showing 
a significant reduction in a number of contribut- 
ing subprocesses and diagrams. 

The method is based on two approximations: 

• mixing with the third generation is ne- 
glected 

VCKM =^ ( 1 ) ' 

y _ f cosdc sindc \ 
y — sin i!)c cos / 



where i?c is the Cabbibo angle. 

• masses of the first four quarks are neglected 
M„ = Md = M,=M, = 

One should stress that usually both approxi- 
mations work very well for collider processes, al- 
though in the case when masses are needed to 
be taken into account one should use the regular 
CompHEP Standard Model. 

Under the assumptions above the method is 
based on the two main ideas: 

• rotation of down quarks, thus, transporting 
the mixing matrix elements from the ver- 
tices of subprocess Feynman diagrams to 
the parton distribution functions 

• the contributing diagrams are splitted into 
exact gauge invariant subclasses |23 for 
each of the subprocess according to a dif- 
ferent topology of the quark lines involved. 

As a result one finds new rules for a convolution 
with the structure functions where each gauge in- 
variant subclass of squared diagrams with certain 
topology of the quark lines is convoluted with 
a corresponding combination of structure func- 
tions. In general there are only three different 
possibilities which we refer to as a convolution 
rules number 1,2 and 3. 

Scattering topology. 1st Rule 

For this topology there are two quark loops af- 
ter squaring, and each loop includes only one ini- 
tial quark. In such a case the structure functions 
are summed over the two flavors for each initial 
state and the dependence on Cabbibo angle dis- 
appears. The following formula gives a particu- 
lar example for this case when one initial state 
is down quarks while another initial state is up 
antiquarks: 

PscP =^ j dxidx2 [fd{xi) + fs{xi)] 
[h{x2) + h{x2)] |A^P 

Annihilation CC topology. 2nd Rule 

For this topology there is only one quark loop 
in the squared diagram which includes both ini- 
tial quarks. For the CC case one of the initial 
state partons is always the up type quark or an- 
tiquark while the another initial state parton is 
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Figure 2. Scattering Topology (see the Figures 
for other topologies in j^) 



the down type quark or antiquark. The effec- 
tive function is the sum of the individual quark 
and antiquark structure function products times 
cos^ -de [siv? -de) in case of the same (different) 
generation for the quark and antiquark. Some 
generic example for the convolution in this case 
is 

=^ j dxidx2[fd{xi)fu{x2)cos''d, 

+fs{xi)fc{x2)cOS^ -de 

+fdixi) Mx2)sin^ 
+fs{xi)fu{x2)sm^^,] 

Annihilation NC topology. 3d Rule 

This topology is similar to the previous one 
leading to one quark loop in the squared diagram. 
But now both initial quark and antiquark are ei- 
ther up or down type. Generic example for the 
convolution with structure functions in this case 
is the following (note, that there is no dependence 
on the Cabbibo angle as for the first rule): 

iPf^P =^ J dxidX2 

[fd{xi)fd{x2)+fs{x^)f-s{x2)] lA^P 

A special rule is used in the cases when the final 
quarks and antiquarks form a loop in a squared 
diagram not involving the initial quarks. In these 
cases each loop of this sort gives a factor of 2 af- 
ter summation over flavors. Of course, this rule 
is valid only under the assumption that the frag- 
mentation of the four light quarks and antiquarks 
leads to indistinguishable jets. If one includes 
nontrivial fragmentation functions, e.g. for a c- 
quark, the above rules have to be modifled. 
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All the above ideas are implemented in Com- 
pHEP in the special model denominated as 
SMuii. For the processes where only QCD inter- 
actions are involved one can simplify further the 
light quark flavor combinatorics since QCD does 
not distinguish different flavors and electroweak 
charges. The corresponding CompHEP model is 
denominated as SMqQ. 

2.2. EVENT COLOR FLOWS 

As will be described below, the CompHEP has 
an interface to PYTHIA program. (The interface 
to HERWIG is in progress). For a fragmentation 
and hadronization PYTHIA needs an information 
on color configurations for the generated partonic 
level events. For the color flows PYTHIA re- 
quires an indication of colorless pairs or, in other 
words, the definition of color chains. In terms of 
QCD the sequences of color chains correspond to 
the elements of the orthogonal color basis in the 
limit of infinite number of colors, Nc — oo. How- 
ever CompHEP, as well as other matrix element 
programs, calculates exactly the matrix element 
squared with the number of colors Nc = 3. 

The procedure described in this section is used 
in CompHEP in order to match the |Mp calcu- 
lation with Nc — 3 and the color chains genera- 
tion needed for further hadronization. Note, that 
originally this algorithm has been implemented in 
CompHEP 4.1.10 by A. Pukhov. 

The matrix element is a sum of Feynman dia- 
grams 

M = J2^dTd 

d 

where is the Lorentz part of the diagram, and 
Td is its color structure . In the limit Nc — oo the 
color structure of diagrams can be decomposed 
over the elements of color basis (color chains) as 

i 

where the elements of the basis are orthonormal, 
U ® tj — Sij. Thus, the squared matrix element 
can be represented in the limit Nc = oo as a sum 
of positive contributions 

\M\%^^^ = J2k^, h = \J2LdCd\^. 

i d 



As a result, for each event (phase space point) 
one gets the set of color chains U with weights 



So the event generation procedure includes two 
steps. First step is the standard generation of an 
event as a phase space point with |Mp calculated 
at = 3 like in the standard event weighting 
procedure. At the second step color chains ti are 
selected with the weighting factors ki, Nc = oo. 

2.3. LANHEP AND SUSY MODELS 

The LanHEP program 25 was developed in 
the framework of the CompHEP project to imple- 
ment easily "new physics" models and/or anoma- 
lous effective Lagrangian terms. LanHEP input is 
the Lagrangian written in the compact form close 
to one used in textbooks. In particular, the La- 
grangian terms as the LanHEP input can be writ- 
ten with summation over indices of the space-time 
or gauge symmetries and using special symbols 
for complicated expressions, such as the covariant 
derivative and the strength tensor for gauge fields. 
There are 2-component spinors and the superpo- 
tential formalism available. Correctness of the 
model can be checked for hermiticity, mass ma- 
trix diagonalization, and BRST invariance. The 
output of the program is a set of Feynman rules 
in terms of physical fields and independent pa- 
rameters. This output can be written in I^TJrjX 
format and in the form of CompHEP model files, 
which allows one to start calculations of processes 
in the new physical model. 

We have used the description of the MSSM 
physical spectrum and the Lagrangian for 3 gen- 
erations of quarks and leptons given in [21]. We 
also add the effective Higgs potential to the tree 
level MSSM Lagrangian to incorporate the radia- 
tive corrections as it is given in the FeynHiggs 
program J57| to the masses and couplings of Higgs 
bosons in a gauge invariant way. Details of the 
notation for the particles and parameters of the 
MSSM in CompHEP are given in [SS] The Feyn- 
man rules generated by LanHEP were success- 
fully compared with the ones listed in as 
well as with many results of other publications. 
The comparison to GRACE/SUSY 2 program 
was successfully done for many processes. 
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The SUGRA and GMSB SUSY models are also 
implemented to the CompHEP by means of link- 
ing with the ISAJET 29 library. There are also 
MSSM extensions available, such as R-parity vi- 
olation model and the model with gravitino and 
sgoldstino fields PH] , 

Lagrangians for other nonstandard models, e.g. 
models with leptoquarks ^31.32^. complete two- 
Higgs-doublet model with conserved or broken 
CP invariance |33| . models with anomalous top 
quark interactions [SJ are also available for Com- 
pHEP calculations. 

2.4. COMPHEP-PYTHIA INTERFACE 

PYTHIA has the built-in database of matrix el- 
ements for hard subprocesses which are basically 
of the 2 2 type. The CompHEP-PYTHIA 
interface (the old version of the interface is de- 
scribed in allows to use the processes 2 
3, 4, 5, 6 computed by means of CompHEP as the 
"new" processes for PYTHIA and to include in 
this way ISR/FSR, hadronization, and decays like 
it is done in PYTHIA. 

CompHEP generates unweighted events and 
writes them to the file events_N.txt where N 
is the number of working session. The com- 
mand mixPEV {syntax: ./mixPEV events-Ni.txt 
eventsJSf2.txt ...) mixes randomly several event 
files corresponding to some subprocesses in one 
event fiow according to their relative weights, 
Oi/TijUj, where Oi is the cross section of the 
i-ih subprocess, and writes events to the file 
Mixed. PEV . This file is ready for input to 
PYTHIA 6.2 and contains the data necessary for 
Les Houches Accord I 15 . 

We provide the interface library libinter- 
face62.a and an example of the mam./ program 
which is very easy to be incorporated to any stan- 
dard PYTHIA main.f program. Using the stan- 
dard Unix "make" command a user links main.o, 
libinterface62.a and pythia62.o to the executable 
file generator.exe . The generator takes events 
from the file Mixed. PEV. The number of re- 
quested events to be generated after PYTHIA is 
defined by user in the file INPARM.DAT 

There are already many examples of the Com- 
pHEP and CompHEP-PYTHIA interface use for 
the real analysis at the Tevatron, LHC and LC. 



Of course, as pointed out in the introduction, 
one has to be careful combining new higher order 
processes calculated with CompHEP-PYTHIA 
interface, and the standard PYTHIA processes. 
Double counting of contributions should be care- 
fully analysed, not only when ISR/FSR are 
switched on, but in many other cases. 

2.5. COLLIDING BEAMS IN COMPHEP 

New option is implemented in CompHEP 4.4 - 
a possibility to enter a beam as the initial state. 
In CompHEP the beam is a set of partons. This 
is a natural terminology for hadronic collisions 
(proton or antiproton beams). In CompHEP 4.4 
this option allows a user to enter also an elec- 
tron, photon or quark itself as a beam. In this 
way one can set, for example, the photon as a 
"parton" in electron (effective photon approxima- 
tion), quarks and gluons as "partons" in a pho- 
ton (resolved photons), W-boson as a "parton" 
in electron or quark (effective W-approximation) 
etc. This option is specially useful for lepton and 
lepton-hadron colliders. Moreover, one can in- 
troduce a beam with arbitrary content for any 
private purposes. 

The parton distribution functions are assigned 
now to the beam by a user before the process is 
entered. All predefined beams are collected into 
a special table which can be modified by a user 
if necessary. When the user sets a process (s)he 
indicates a particle or a beam (from the table) as 
an initial state. The corresponding table can be 
displayed on the screen by pressing the F3 key. 

In the CompHEP 4.4 the latest versions of the 
CTEQ parton distribution functions |35j are im- 
plemented. However, some previous sets are also 
available since for an analysis where they could 
be needed. 

2.6. SYMBOLIC BATCH MODE 

The main goal of the symbolic batch script is 
to launch the symbolic calculation of squared di- 
agrams using a command line interface (the non- 
GUI mode). This allows the reuse of the script, 
avoiding typing mistakes, and also helps to send 
large symbolic calculation tasks to the computer 
farms. 

The script imitates the typing input data in 
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CompHEP menu (initial states, energies, final 
particles etc) . This option was described in ^| - 
so-called blind mode in CompHEP 4.1.10. 

All the data, necessary for further symbolic cal- 
culations, are taken from a file with a default 
name process.dat. Detailed instructions and ex- 
planations for the parameters are contained in 
this file. 

After filling in the file process.dat a user 
launches the script. The results of the script ex- 
ecution will be the binary file s_comphep.exe and 
the file symb_batch.log where all details on the bi- 
nary creation are stored. 

One can set another name for the data file or 
the "results" directory by the -n and -d options 
respectively. If the CompHEP directory " results" 
is not empty before the script is launched, it is 
renamed to results_old_0 and the script creates 
a new "results" directory. The symbolic batch 
script has some extra options: 

-recovery if the s_comphep program 
crashes during calculations of squared di- 
agrams, one can launch the script with 
the recovery option, and the s_comphep re- 
sumes computations from the last calcu- 
lated diagram; 

-relink if the user has changed the 
userFun.c file and wants to relink the 
n_comphep program with a new file, (s)he 
can launch the script with the relink option. 
All details of the relinking are stored in the 
symb_batch.log file; 

-show diagrams this option is applied if 
the user wants to exclude by hand some di- 
agrams from the calculation. The script 
launches s_comphep in GUI mode in the 
Feynman diagrams menu. After reviewing 
the diagrams the user has to finish the GUI 
session and the script will go on. 

2.7. NUMERICAL BATCH MODE 

The CompHEP was originally created as a pro- 
gram with an interactive GUI interface. However 
in practice for many cases a command line in- 
terface (non-GUI) is very useful and it was im- 
plemented in the CompHEP 4.4 numerical batch 



mode. The interface allows to perform large time 
scale calculations using, if necessary, the com- 
puter farms (PBS and LSF systems), to compute 
processes with a large number of subprocesses, in 
particular, in parallel if computer resources are 
enough. The interface is very useful for repeating 
numerical calculations by non-expert users with- 
out making trivial mistakes. 

In this mode numerical computations start af- 
ter finishing the symbolic computations and cre- 
ation of the executable file n_comphep.exe in the 
directory "results". At the first step a user has 
to create the numerical batch data file batch.dat 
by launching the command . / numjjatch.pl from 
the user working directory, where CompHEP has 
been started. The created file batch.dat contains 
the default parameters (model parameters, cuts, 
kinematics, etc.) listed in the file session.dat in 
the directory "results". If necessary the user can 
change any of these parameters by editing the 
file batch.dat or using GUI interface to modify 
initial session.dat file. In the latter case one can 
start the command ./num -batch. pi with the op- 
tion add-ses2bat to include parameters from 

the new session.dat file to the batch.dat. 

The command (script) ./numjjatch.pl has 
many other useful options, for example: 

—help - to list the complete set of options 
with an explanation 

-d dir_name - to change the default direc- 
tory results to the directory dir_name 
— add_ses2bat [file] - to add parameters 
from existing file session.dat in the direc- 
tory results or from [file], if indicated, to 
the batch.dat file 

-run [vegas,max,evnt] - to start numer- 
ical calculations including all steps by de- 
fault or only the steps as indicated by the 
keys [vegas,max,evnt], where [vegas] means 
the cross section calculation, [max] - search 
for maxima in each of the phase space sub- 
cubes, [evnt] - unweighted event generation 

-run [cleanstat,cleangrid, clean] - to 

clear statistics or grid or both 

-proc nl,n2-n3,n4... - to run only the 
subprocesses with corresponding numbers 
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-pbs ["pbs prefix"] - to start calculations 
in parallel with PBS batch system. The de- 
fault PBS prefix is: " qsub -I " 
-Isf ["Isf prefix"] - to start calculations in 
parallel with LSF batch system. The de- 
fault LSF prefix is: "bsub -I " 

For example, the command 
./ numJbatch.pl — run vegas —pbs 
creates temporary subdirectories for each of the 
subprocess and run the PBS batch job for each of 
subprocess with PBS command: qsub -I run.sh, 
where the script run.sh is created automatically 
and starts the n_comphep.exe. These temporary 
subdirectories will be automatically removed and 
the results for all subprocess will be saved in one 
directory. 

A user can also calculate subprocesses in par- 
allel without the PBS (LSF) batch system by 
starting the same command with empty PBS 
(LSF) prefix: 

./num_batch.pl -run vegas -Isf ' ' 

2.8. FUTURE PLANS 

There are several directions for the future de- 
velopment of CompHEP. One important direc- 
tion is a further development of distributed Monte 
Carlo calculation and event generation on com- 
puter clusters |H7j . Another important direc- 
tion is an implementation of the FORM com- 
puter algebra program f3H1 for symbolic calcula- 
tions of matrix elements (some preliminary real- 
ization was described in f3Sj). This option will 
allow to introduce new complicated structures in 
the vertices (e.g. form- factors), implement new 
algorithms to increase the efficiency of symbolic 
calculations, perform calculations in theories with 
extra dimensions, use the dimensional rcgulariza- 
tion, perform polarized calculations by introduc- 
ing the corresponding density matrices for exter- 
nal lines of squared diagrams. We plan to cre- 
ate the Les Houches accord 1 based interface to 
HERWIG and the Les Houches accord 2 based in- 
terface to PDFs as well as the SUSY Les Houches 
accord based interface 001 to various codes cal- 
culating the SUSY mass spectrum. 

The long term CompHEP collaboration 
projects which are under discussion include sym- 



bolic and numerical amplitude calculations with 
extension to the 1-loop case using various meth- 
ods, incorporation the gauge invariant classes of 
diagrams, etc. 

3. THE COMPHEP BASED GENERA- 
TOR SINGLETOP 

The CompHEP has been used to prepare a spe- 
cial event generator SingleTop to simulate the 
electroweak single top quark production with its 
subsequent decays at the Tevatron and LHC. Sin- 
gle top is expected to be discovered at the Teva- 
tron Run II and will a very interesting subject 
of detail studies at the LHC (see the reviews 

mm)- 

There are three main processes of singe top 
production at hadron colliders characterizing by 
the virtuality of the participating at the pro- 
cess ly-boson: t-channel, s-channel and associ- 
ated tW mechanisms respectively. 

The generator SingleTop includes all the three 
processes and provides Monte-Carlo unweighted 
events at the NLO QCD level. We discuss shortly 
here only the main process with the largest rate, 
the t-channel production. The representative 
NLO diagrams are shown in Fig|31 The top de- 
cay is not shown, however it is included with all 
the spin correlations. 




1.1 1.2 1.3 1.4 



Figure 3. LO order and representative loop and 
tree NLO diagrams to the t-channel single top 
production 



We compute by means of the CompHEP the 
LO order process 2^2 with the b-quark in the 
initial state and top spin correlated 1 — > 3 sub- 
sequent decay, put it into PYTHIA using the in- 
terface and switch on ISR/FSR. Then with Com- 
pHEP we compute the NLO tree level corrections 
-2^3 processes with additional b- and light 
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quarks or gluons in the final state including also 
the top decay with spin correlations. We split 
the phase space region in " soft" and " hard" parts 
on pt of those additional b and light jets being 
from PYTHIA radiation in the "soft" and from 
the CompHEP matrix element calculation in the 
"hard" regions. The soft part is normalized in 
such a way that all parts being taken together give 
known from calculations the NLO cross section 
|43I44| . The splitting parameters are turned such 
that all the distribution become smooth after the 
normalization. The performed cross checks show 
that the computed NLO distributions 01] are cor- 
rectly reproduced. Therefore, prepared in that 
way generator does not have a double counting, 
gives the NLO rate and distributions, and include 
all the spin correlations. 

The first release of the generator jlHl did not 
include the hard radiation of the light jets, while 
the latest version [3^ currently used in the anal- 
ysis by the Fermilab DO and the LHC CMS col- 
laborations includes all the mentioned properties. 

There are several recent examples of the Com- 
pHEP use for various simulations, namely, the 
generator for MSSM Higgs bosons in the intense 
coupling regime at the LHC 07] , the generator for 
sfermion pair production with their subsequent 
decays to polarized fermions ^IS^ and the genera- 
tor for measurements of the HZj coupling of jlH] 
at Linear Colliders. 

4. CONCLUSIONS 

The CompHEP package is a powerful tool for 
a simulation of different physical processes at 
hadron and lepton colliders. The most important 
CompHEP 4.4 improvements include the build- 
in MSSM, simplification of quark flavor combi- 
natorics, generation of unweighted events with 
color chains defined for the following showering 
and hadronization procedure, Les Houches ac- 
cord based interface with PYTHIA, symbolic and 
numerical batch modes, the latest PDF sets im- 
plementation and a new treatment of the initial 
beams. 

The CompHEP is basically the LO program. 
However it allows to include partly the NLO cor- 
rections. NLO tree level 2 ^ -I- 1 correc- 



tions to the 2 ^ N process can be computed. 
One can include the NLO structure functions and 
loop relations between the parameters like in the 
MSSM, K-factors when available, and loop con- 
tributions from the existing publications as form- 
factors. However a correct matching to complete 
NLO is obviously a nontrivial problem which has 
to be carefully considered in each particular phys- 
ical case. The CompHEP based generator Single- 
Top is an example of a reasonable matching of 
the different-order contributions. 
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